The unusual gene architecture of polyubiquitin is created by dual-specific splice sites.

BACKGROUND: The removal of introns occurs through the splicing of a 5' splice site (5'ss) with a 3' splice site (3'ss). These two elements are recognized by distinct components of the spliceosome. However, introns in higher eukaryotes contain many matches to the 5' and 3' splice-site motifs that are presumed not to be used. RESULTS: Here, we find that many of these sites can be used. We also find occurrences of the AGGT motif that can function as either a 5'ss or a 3'ss-previously referred to as dual-specific splice sites (DSSs)-within introns. Analysis of the Sequence Read Archive reveals a 3.1-fold enrichment of DSSs relative to expectation, implying synergy between the ability to function as a 5'ss and 3'ss. Despite this suggested mechanistic advantage, DSSs are 2.7- and 4.7-fold underrepresented in annotated 5' and 3' splice sites. A curious exception is the polyubiquitin gene UBC, which contains a tandem array of DSSs that precisely delimit the boundary of each ubiquitin monomer. The resulting isoforms splice stochastically to include a variable number of ubiquitin monomers. We found no evidence of tissue-specific or feedback regulation but note the 8.4-fold enrichment of DSS-spliced introns in tandem repeat genes suggests a driving role in the evolution of genes like UBC. CONCLUSIONS: We find an excess of unannotated splice sites and the utilization of DSSs in tandem repeats supports the role of splicing in gene evolution. These findings enhance our understanding of the diverse and complex nature of the splicing process.

The debranching enzyme Dbr1 regulates lariat turnover and intron splicing.

The majority of genic transcription is intronic. Introns are removed by splicing as branched lariat RNAs which require rapid recycling. The branch site is recognized during splicing catalysis and later debranched by Dbr1 in the rate-limiting step of lariat turnover. Through generation of the first viable DBR1 knockout cell line, we find the predominantly nuclear Dbr1 enzyme to encode the sole debranching activity in human cells. Dbr1 preferentially debranches substrates that contain canonical U2 binding motifs, suggesting that branchsites discovered through sequencing do not necessarily represent those favored by the spliceosome. We find that Dbr1 also exhibits specificity for particular 5' splice site sequences. We identify Dbr1 interactors through co-immunoprecipitation mass spectroscopy. We present a mechanistic model for Dbr1 recruitment to the branchpoint through the intron-binding protein AQR. In addition to a 20-fold increase in lariats, Dbr1 depletion increases exon skipping. Using ADAR fusions to timestamp lariats, we demonstrate a defect in spliceosome recycling. In the absence of Dbr1, spliceosomal components remain associated with the lariat for a longer period of time. As splicing is co-transcriptional, slower recycling increases the likelihood that downstream exons will be available for exon skipping.

Activation of human RNA lariat debranching enzyme Dbr1 by binding protein TTDN1 occurs though an intrinsically disordered C-terminal domain.

In eukaryotic cells, the introns are excised from pre-mRNA by the spliceosome. These introns typically have a lariat configuration due to the 2'-5' phosphodiester bond between an internal branched residue and the 5' terminus of the RNA. The only enzyme known to selectively hydrolyze the 2'-5' linkage of these lariats is the RNA lariat debranching enzyme Dbr1. In humans, Dbr1 is involved in processes such as class-switch recombination of immunoglobulin genes, and its dysfunction is implicated in viral encephalitis, HIV, ALS, and cancer. However, mechanistic details of precisely how Dbr1 affects these processes are missing. Here we show that human Dbr1 contains a disordered C-terminal domain through sequence analysis and nuclear magnetic resonance. This domain stabilizes Dbr1 in vitro by reducing aggregation but is dispensable for debranching activity. We establish that Dbr1 requires Fe2+ for efficient catalysis and demonstrate that the noncatalytic protein Drn1 and the uncharacterized protein trichothiodystrophy nonphotosensitive 1 directly bind to Dbr1. We demonstrate addition of trichothiodystrophy nonphotosensitive 1 to in vitro debranching reactions increases the catalytic efficiency of human Dbr1 19-fold but has no effect on the activity of Dbr1 from the amoeba Entamoeba histolytica, which lacks a disordered C-terminal domain. Finally, we systematically examine how the identity of the branchpoint nucleotide affects debranching rates. These findings describe new aspects of Dbr1 function in humans and further clarify how Dbr1 contributes to human health and disease.

A functional link between lariat debranching enzyme and the intron-binding complex is defective in non-photosensitive trichothiodystrophy.

The pre-mRNA life cycle requires intron processing; yet, how intron-processing defects influence splicing and gene expression is unclear. Here, we find that TTDN1/MPLKIP, which is encoded by a gene implicated in non-photosensitive trichothiodystrophy (NP-TTD), functionally links intron lariat processing to spliceosomal function. The conserved TTDN1 C-terminal region directly binds lariat debranching enzyme DBR1, whereas its N-terminal intrinsically disordered region (IDR) binds the intron-binding complex (IBC). TTDN1 loss, or a mutated IDR, causes significant intron lariat accumulation, as well as splicing and gene expression defects, mirroring phenotypes observed in NP-TTD patient cells. A Ttdn1-deficient mouse model recapitulates intron-processing defects and certain neurodevelopmental phenotypes seen in NP-TTD. Fusing DBR1 to the TTDN1 IDR is sufficient to recruit DBR1 to the IBC and circumvents the functional requirement for TTDN1. Collectively, our findings link RNA lariat processing with splicing outcomes by revealing the molecular function of TTDN1.

Metal content and kinetic properties of yeast RNA lariat debranching enzyme Dbr1.

In eukaryotic cells, intron lariats produced by the spliceosome contain a 2'5' phosphodiester linkage. The RNA lariat debranching enzyme, Dbr1, is the only enzyme known to hydrolyze this bond. Dbr1 is a member of the metallophosphoesterase (MPE) family of enzymes, and recent X-ray crystal structures and biochemistry data demonstrate that Dbr1 from Entamoeba histolytica uses combinations of Mn2+, Zn2+, and Fe2+ as enzymatic cofactors. Here, we examine the kinetic properties and metal dependence of the Dbr1 homolog from Saccharomyces cerevisiae (yDbr1). Elemental analysis measured stoichiometric quantities of Fe and Zn in yDbr1 purified following heterologous expression E. coli We analyzed the ability of Fe2+, Zn2+, and Mn2+ to reconstitute activity in metal-free apoenzyme. Purified yDbr1 was highly active, turning over substrate at 5.6 sec-1, and apo-yDbr1 reconstituted with Fe2+ was the most active species, turning over at 9.2 sec-1 We treated human lymphoblastoid cells with the iron-chelator deferoxamine and measured a twofold increase in cellular lariats. These data suggest that Fe is an important biological cofactor for Dbr1 enzymes.

Massively parallel reporter assays discover de novo exonic splicing mutants in paralogs of Autism genes.

To determine the contribution of defective splicing in Autism Spectrum Disorders (ASD), the most common neurodevelopmental disorder, a high throughput Massively Parallel Splicing Assay (MaPSY) was employed and identified 42 exonic splicing mutants out of 725 coding de novo variants discovered in the sequencing of ASD families. A redesign of the minigene constructs in MaPSY revealed that upstream exons with strong 5' splice sites increase the magnitude of skipping phenotypes observed in downstream exons. Select hits were validated by RT-PCR and amplicon sequencing in patient cell lines. Exonic splicing mutants were enriched in probands relative to unaffected siblings -especially synonymous variants (7.5% vs 3.5%, respectively). Of the 26 genes disrupted by exonic splicing mutations, 6 were in known ASD genes and 3 were in paralogs of known ASD genes. Of particular interest was a synonymous variant in TNRC6C - an ASD gene paralog with interactions with other ASD genes. Clinical records of 3 ASD patients with TNRC6C variant revealed respiratory issues consistent with phenotypes observed in TNRC6 depleted mice. Overall, this study highlights the need for splicing analysis in determining variant pathogenicity, especially as it relates to ASD.

Efficient Detection of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) from Exhaled Breath.

Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) is transmitted through airborne particles in exhaled breath, causing severe respiratory disease, coronavirus disease-2019 (COVID-19), in some patients. Samples for SARS-CoV-2 testing are typically collected by nasopharyngeal swab, with the virus detected by PCR; however, patients can test positive for 3 months after infection. Without the capacity to assay SARS-CoV-2 in breath, it is not possible to understand the risk for transmission from infected individuals. To detect virus in breath, the Bubbler-a breathalyzer that reverse-transcribes RNA from SARS-CoV-2 particles into a sample-specific barcoded cDNA-was developed. In a study of 70 hospitalized patients, the Bubbler was both more predictive of lower respiratory tract involvement (abnormal chest X-ray) and less invasive than alternatives. Samples tested using the Bubbler were threefold more enriched for SARS-CoV-2 RNA than were samples from tongue swabs, implying that virus particles were being directly sampled. The barcode-enabled Bubbler was used for simultaneous diagnosis in large batches of pooled samples at a lower limit of detection of 334 genomic copies per sample. Diagnosis by sequencing can provide additional information, such as viral load and strain identity. The Bubbler was configured to sample nucleic acids in water droplets circulating in air, demonstrating its potential in environmental monitoring and the protective effect of adequate ventilation.

Mutational bias and the protein code shape the evolution of splicing enhancers.

Exonic splicing enhancers (ESEs) are enriched in exons relative to introns and bind splicing activators. This study considers a fundamental question of co-evolution: How did ESE motifs become enriched in exons prior to the evolution of ESE recognition? We hypothesize that the high exon to intron motif ratios necessary for ESE function were created by mutational bias coupled with purifying selection on the protein code. These two forces retain certain coding motifs in exons while passively depleting them from introns. Through the use of simulations, genomic analyses, and high throughput splicing assays, we confirm the key predictions of this hypothesis, including an overlap between protein and splicing information in ESEs. We discuss the implications of mutational bias as an evolutionary driver in other cis-regulatory systems.